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1. Introduction 

This paper is based on a lecture given at the workshop “Flavorful Ways to New Physics". However, 
its scope extends beyond the topics that could be covered in the limited time of the lecture. The aim 
is to give a broad but at the same time sufficiently in-depth overview of the topic of charm physics 
to be of benefit to students with a minimum of background in particle physics. The structure 
and content of the paper is based on a review article [1], with several additional topics, more 
pedagogical detail as well as a content that has been brought up to date. 

Charm physics covers the studies of a range of composite particles containing charm quarks 
which provide unique opportunities for probing the strong and weak interactions in the standard 
model and beyond. The charm quark, being the up-type quark of the second of the three gener¬ 
ations, is the third-heaviest of the six quarks. Charm particles can exist as so-called open charm 
mesons or baryons, containing one or several (for baryons) charm quarks, or as charmonium states 
which are bound states of charm anti-charm (cc) quark pairs. In addition, several states containing 
charm quarks have been observed that cannot be explained as conventional mesons or baryons or 
as bound states thereof. Their exact nature is one of today’s unanswered questions. 

The uniqueness of charm particles lies in their decays. The charm quark can only decay via 
weak decays, mediated by a W ± -boson, into a strange or down quark. An exception to this are 
decays of ground state charmonium mesons, which decay via annihilation of the charm and anti¬ 
charm quarks. Thus, open charm particles are the only ones allowing the study of weak decays of 
an up-type quark in a bound state. 

In 2009, Ikaros Bigi asked whether charm’s third time could be the real charm [2]. Charm’s 
first time was the discovery of the J/iy [3, 4], which followed three years after the possible first 
observation of an open charm decay in cosmic ray showers [5]. This discovery confirmed the 
existence of a fourth quark as expected by the GIM mechanism [6] motivated by the non-existence 
of flavour-changing neutral currents [7, 8] in conjunction with the observation of the mixing of 
neutral kaons [9, 10, 11]. The second time charm attracted considerable attention was caused by 
the observation of D s j states, [12, 13, 14, 15] which could not be accommodated by QCD motivated 
quark models [16, 17, 18, 19, 20]. Until today, excited charmonium and open charm particles 
provide an excellent laboratory for studying QCD. 

Charm’s third time started with the first evidence for mixing of neutral charm mesons reported 
by BaBar [21] and Belle [22] in 2007. Since then a lot of work went into more precise mea¬ 
surements of the mixing phenomenon as well as into searches for charge-parity (CP) symmetry 
violation in the charm sector. At the same time theoretical calculations were improved even though 
precise standard model predictions are still a major challenge. 

In 2013, the observation of a charged charmonium-like state has marked the beginning of 
what might be called charm's fourth time. The state with a mass of about 3900MeV/c 2 and a 
width in the region of 50MeV/c 2 was first observed in decays to J/\fjn + final states [23, 24]. A 
similar state or possibly the same has been observed at the threshold of the DD* spectrum [25]. Its 
charmonium content and electrical charge rule out any interpretation as a two-quark or three-quark 
state. However, also attempts to identify the state e.g. as a molecule of two charm mesons or as 
a tetraquark state do not yield clear results. Already in 2008, the Belle collaboration observed a 
similar charged charmonium-like state around a mass of 4430MeV/c 2 [26]. However, this could 
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not be confirmed by a BaBar analysis [27]. Further studies by the Belle collaboration [28, 29] gave 
additional insight into spin and parity of the state and supported their initial measurement. This 
was finally confirmed by an LHCb analysis based on a significantly larger dataset [30]. Detailed 
reviews of these states are available in Refs. [31, 32]. 

The focus of this paper lies on the current situation of studies of processes, which are mediated 
by the weak interaction, using open charm particles. Particular focus is given to mixing and CP 
violation, followed by comments on rare charm decays at the end of this review. 

1.1 Charm production 

Charm physics has been and is being performed at a range of different accelerators. These come 
with different production mechanisms and thus with largely varying production cross-sections. At 
e 1 e colliders two different running conditions are of interest to charm physics. Tuning the centre - 
of-mass energy to slightly more than 3770MeV/c 2 resonantly produces t/r(3770) mesons. These 
decay almost exclusively to quantum-correlated D°D° or D + D~ pairs. The use of this produc¬ 
tion type has been pioneered by the MARK III collaboration at the e 1 e ~ storage ring SPEAR at 
SLAC [33]. More recently, this production mode has been used for the CLEO-c experiment at the 
CESR-c collider as well as for BESIII at BEPCII. For ah these experiments the collisions happen 
at rest, which removes the possibility to study the decay-time structure of the D mesons. 

The most commonly used alternative is running at a higher centre-of-mass energy to produce 
T(4S) on resonance which decay into quantum-correlated B°B° or B B pairs. This is used by the 
BaBar and Belle experiments which are located at the PEP-II and KERB colliders, respectively. 
Both PEP-II and KERB are asymmetric colliders, which means that the energies of the e + and e 
beams differ. Flence, they have a collision system that is boosted with respect to the laboratory 
frame. This allows measurements with decay-time resolutions about a factor two to four below 
the D° lifetime and therefore decay-time dependent studies. Currently, the Super-KERB collider is 
under construction, which will act as a 71-factory with much improved luminosity for the Belle II 
experiment. 

The production cross-section for producing DD pairs at the i/r(3770) resonance is approxi¬ 
mately 8nb [34], When running at the Y(4S) resonance, the cross-section for producing at least 
one D° meson is 1.45nb [35]. The latter scenario gives access to all species of charm particles 
while the i//(3770) only decays into D°D° or D + D~ pairs. The BaBar and Belle experiments have 
collected integrated luminosities of about 500 lb 1 and 1000 lb 1 , respectively. The Belle II exper¬ 
iment is expected to collect 50ab _1 between 2017 and 2022. CLEO-c has collected 0.5 fb _1 at the 
t/r(3770) resonance as well as around 0.3 fb 1 above the threshold for 7) s ' 7) s production. BESIII 
has so far collected nearly 3fb~' at the i/r(3770) resonance as well as nearly 6fb _1 at energies 
between 3.9GeV and 4.6GeV. 

At hadron colliders the production cross-sections are significantly higher. The cross-section 
for producing cc pairs in proton-proton collisions at the LHC with a centre-of-mass energy 7TeV 
is about 6nrb [36], i.e. more than six orders of magnitude higher compared to operating an e 1 e 
collider at the Y(4S) resonance. This corresponds to a cross-section of about 1.4 mb for producing 
D° in the LHCb acceptance 1 [37], This number may be compared to its equivalent at CDF which 

'Given as a range in momentum transverse to the beam direction and rapidity as pj < 8GeV/c,2 < y < 4.5. 
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Experiment 

Year 

vA 

Oacc{D°) 

L 

n(D°) 

CLEO-c 

2003-2008 

3.77 GeV 

8nb 

0.5 flG 1 

4.0 x 10 6 

BESIII 

2010-2011 

3.77 GeV 

8nb 

3 fb 1 

2.4 x 10 7 

BaBar 

1999-2008 

10.6 GeV 

1.45 nb 

500fb" 1 

7.3 x 10 8 

Belle 

1999-2010 

10.6-10.9 GeV 

1.45 nb 

1000fb- 1 

1.5 x 10 9 

CDF 

2001-2011 

2 TeV 

13 /Tb 

10ft)- 1 

1.3 x 10 11 

LHCb 

2011 

7 TeV 

1.4 mb 

1 fb 1 

1.4 x 10 12 

LHCb 

2012 

8 TeV 

1.6mb* 

2 fb 1 

3.2 x 10 12 


Table 1: Charm production values for different experiments based on the example of the production of D° 
mesons inside the respective detector acceptances as defined in the text. The LHCb D° production cross- 
section at y/s = 8 TeV has been extrapolated from that at y/s = 7 TeV assuming linear scaling with y/s. 

has been measured to 13 /f b inside the detector acceptance 2 for proton anti-proton collisions at the 
Tevatron with y/s = 1.96TeV [38]. CDF has collected a total of about lOflG 1 while LHCb has 
collected 3 fb 1 during the first run of the LHC, corresponding to0.13xl0 12 and 4.6 x 10 12 D° 
mesons produced in the respective detector acceptances. Table 1 summarises the different charm 
productions with the example of the production of D° mesons. 

The production of charm quarks in hadron collisions occurs predominantly in very asymmet¬ 
ric collisions which result in heavily boosted quarks with high rapidities. While the asymmetry of 
the collisions has its origin in two different beam energies at 5-factories, it comes from the fact 
that LHC collisions occur among partons rather than the protons as a whole. At LHC energies it is 
mostly the gluons participating in the collisions and the fractional energy of the two colliding par- 
tons is likely to differ. This difference leads to a smaller centre-of-mass energy as well as a boost 
of the collision system. Therefore, LHCb is ideally suited for performing decay-time dependent 
studies of charm decays. At the same time, the cc cross-section at the LHC is about 10% of the 
total inelastic cross-section which allows to have reasonably low background levels for a hadronic 
environment. The coverage of nearly the full solid angle of the e + e~ -collider experiments men¬ 
tioned here makes them very powerful instruments for analysing decays involving neutral particles 
that may remain undetected or for inclusive studies. In addition, the low multiplicity of collision 
products makes these collisions particularly clean laboratories due to the low level of backgrounds. 

1.2 Charm spectroscopy 

Open charm mesons and baryons and charmonium mesons exhibit a rich structure of excited states. 
Based on the quark model the expected spectra have been predicted thirty years ago [16]. The IS 
ground states D^, and D u have spin-parity J r = 0 , and 1~, respectively, and exist for the neutral 
(cu) and charged (cd, cs) mesons. In general, the configuration with 5 = (— I ) J is called natural 
parity and denoted as D* , while P = — (— I ) J is called unnatural parity. The ground state transition 
D*+ —> D°7i + will be of importance for identifying the flavour of D° mesons later in this paper. 

The first orbital excitations (nL = IP multiplets) have been already observed for both charm 
and charm-strange mesons. For the latter for example, these are D* 0 (2317) (0+), two 1+ states 
D v ] (2460) and D s i (2536), as well as D* 2 (2573) with spin-parity 2+. The two lightest ones were 

2 The CDF acceptance is defined as pj > 5.5GeV/c, [y| < 1. 
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found to be much lighter and narrower than the quark models predicted. This triggered a lively dis¬ 
cussion about the nature of these mesons and led to the revision of the corresponding calculations. 
There are also a number of candidates for the first radial excitations (n = 2) and/or higher orbital 
excitations observed. 

A similarly rich spectrum exists for charmonium, i.e. cc mesons. The ground states are rj c (l S) 
and J/\j/(lS) with spin-parity J p = 0~, and 1~, respectively. A detailed review of charmonium 
spectroscopy can be found in Ref. [39]. Thirdly, charm baryons also have a rich spectroscopy; 
however, this is the least explored sector as so far only singly charmed baryons have been observed. 

A number of the quark model states have now been experimentally confirmed. However, there 
are observed states that cannot be matched to spectroscopy predictions. In addition to spectroscopy, 
the study of production and decay modes of the individual states are needed to pin down their 
precise nature. More exotic models exist that attempt their explanation, such as meson molecules, 
quark-gluon hybrids, glueballs. Other approaches are based particles made of more than three 
quarks to form tetraquarks, pentaquarks, or H-dibaryons. 


1.3 Mixing of neutral mesons 

For neutral mesons, the mass eigenstates, i.e. the physical particles of defined mass and lifetime, 
do not a priori coincide with the flavour eigenstates, which are being produced and evolve in time 
into their anti-particles to finally decay as mass eigenstates. In general the state of a meson can be 
expressed as a linear combination of the flavour eigenstates as 

\ ¥ )=a(t)\M Q )+b(t)\M°), (1.1) 


which is conveniently abbreviated to 


\¥) = 



(1.2) 


For the purpose of mixing it is useful to apply the Weisskopf-Wigner approximation [40,41] by ne¬ 
glecting the flavour-conserving strong and electromagnetic interactions and writing the Schrodinger 
equation 


= Jtfy, 


(1-3) 


with the effective Hamiltonian given by 


j ? = m - -r = 
2 


( M n 
\M 2l 


^Ft 1 M\2 — ^T\2 A 

^21 m 22 - ^ r 22 J ’ 


(1.4) 


where the mass matrix M and the decay matrix F are Hermitian. The eigenstates of the Hamiltonian, 
|Afi, 2 >, have eigenvalues 


Ai, 2 = »F,2 — 2 r i,2> 


(1.5) 


as well as the observable masses m\i and decay widths T] ?■ The time evolution of the physical 
states is therefore given as 


\M {p _(t)) = e~ imi ' 2t e~ Tl ' 2t / 2 \M\ , 2 (0)). 


( 1 . 6 ) 
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Meson 

Mass in GeV 

Width in ps 1 

Lifetime in ps 

X 

y 

K° 

0.49761(2) 

0.005594(2) 

89.54(4) 

51160(210) 

0.946(2) 

0.997(1) 

D° 

1.86484(7) 

2.438(9) 

0.4101(15) 

0.0039(+}^) 

0.0065 (l[g) 

K 

5279.58(17) 

0.658(2) 

1.519(5) 

0.774(6) 

<0.09 


5366.77(24) 

0.869(4) 

1.512(7) 

26.85(13) 

0.069(6) 


Table 2: Overview of parameters relevant to meson mixing. Masses and widths are the average of the two 
physical eigenstates. The uncertainties are given in brackets as multiples of the least significant digit. The 
values used are taken from Ref. [44]. 


Assuming CPT symmetry, the physical eigenstates can be expressed as 


\M u2 )=p\M°)±q\M °>, 


(1.7) 


with complex coefficients satisfying \p\ 2 + \q\ 2 = 1. The choice of the phase convention of p and 
q can be made such that, in the limit of no CP violation, i.e. \q/p\ = 1, CP\M°) = —|M°). This is 
equivalent to stating that \M\) is CP odd and \M 2 ) is CP even. 

Solving Eq. 1.7 for the flavour eigenstates one can express their time dependence by inserting 
Eq. 1.6 and using Eq. 1.7 for the initial states to obtain 


\M°(t)) =/+(f)|M°) + —/_(f)|M°) 

I M°(0) = Mt)\M°) + Pf_(t)\M°), (1.8) 

9 


with 


/±(0 


1 

2 


e~ im 


e -rv / 2 



,—iAmt AYt/l\ 
)' 


(1.9) 


where the quantities Am = m 2 — mi and Ar = Tt — T| have been introduced. It is useful to in¬ 
troduce further the corresponding averages m = ( m\ + m 2 )/2 and T = (T| +r2)/2 as well as the 
dimensionless quantities x = Am/T and y = Ar / (2r). 

At this stage the transition probability of one flavour eigenstate to a state of the same or oppo¬ 
site flavour can be calculated as 


P(M°(t) M°) = P(M°(t) —)• M°) 

P(M°(t) —> M°) 

P(M°(t) —^ M°) 


I/+WI 2 = ^ 


Tt 


\f-{t)\ 2 =\ % 

l/-(0| 2 = 2 f 


(cosh(yr?) + cos(vr?)), 

2 


e 1 '(cosh(yr/‘) — cos(xr?)), (1.10) 


,,-n 


(cosh(yrt) — cos(xrt)). 


These equations show that it is the mixing parameters x, and y, which define the characteristic be¬ 
haviour of neutral meson mixing. In particular it is the mass difference of the physical eigenstates, 
which determines the rate of oscillations between mesons and anti-mesons. Different decay widths 
lead to a non-oscillating change of the exponential decay. A more detailed introduction can be 
found e.g. in Refs. [42, 43]. 

The four systems, which are subject to mixing, are kaons ( K ), charm (D), /?[[. and B ( * mesons. 
For charm mesons the mixing parameters are drastically different compared to those of kaons or 
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Figure 1: The widths and mass differences of the physical states of the flavoured neutral mesons. The 
width corresponds to the inverse lifetime while the mass difference determines the oscillation frequency. 
The values used are taken from Ref. [44]. 


B mesons. Their masses, lifetimes, and mixing parameters are shown in Table 2. Figure 1 gives a 
graphical representation of the widths and mass differences of the four neutral meson systems. The 
kaon system is the only one to have y ~ 1 , resulting in two mass eigenstates with vastly different 
lifetimes, hence their names A'-short (Kg) and A'-long (Aj°). Furthermore, also x ~ I which results 
in a sizeable sinusoidal oscillation frequency as shown in Eq. 1.10). The two A-meson systems 
have reasonably small width splitting; however, they have sizeable values for x. Particularly for the 
B 5 system this leads to fast oscillations which require high experimental accuracy to be resolved. 
Similarly, only recent measurements were able to measure the B® width splitting to be unambigu¬ 
ously non-zero. The charm meson system is the only one where both x and y are significantly less 
than 1, hence the nearly overlapping curves in Fig. 1. 

Experimentally, the different mixing parameters lead to rather different challenges for mea¬ 
surements in the various meson systems. The vast lifetime difference in the kaon system leads to 
the possibility of studying nearly clean samples of just one of the two mass eigenstates by either 
measuring decays close to a production target where K<1 decays dominate, or far away where most 
Kg have decayed before entering the detection region. In the B systems the oscillation frequency 
puts a challenge to the decay-time resolution, particularly for B® mesons as mentioned before. The 
smallness of y requires, to first order, large data samples to acquire the necessary statistical pre- 
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cision for measuring such a small quantity. The latter is particularly true for the charm system, 
where both x and y are small. This is the reason why it was only in 2007 when first evidence for 
charm mixing was observed and took until 2013 for mixing to be observed with high statistical 
significance in a single measurement. 


1.4 CP violation 


The symmetry under CP transformation, i.e. the exchange of particles and anti-particles can be 
violated in different ways. One is a difference in the transition probability of mesons to anti¬ 
mesons compared to the reverse process, i.e. P(M°(t) —> M°) ^ P(M°(t) M°). As can be seen in 

Eq. 1.10, this is the case if \q/p\ / 1. This type of CP violation is independent of the decay mode. 

When considering also the decay of a meson, CP violation is given by the inequality of the 
magnitude of the decay amplitude of a meson to a final state /, |A/|, compared to the magnitude 
of the anti-meson decay amplitude to the CP conjugate final state f, \Aj-\. Hence, CP violation is 
present when \Aj/Af \ / 1. This decay-mode dependent CP violation is called direct CP violation, 
with the corresponding parameter defined as 


_ \ A f\ 2 — \Af\ 2 

“ i= \C?+W 


(Ml) 


Finally, for neutral meson decays the mixing and decay amplitudes can interfere, which can 
lead to a third type of CP violation. This CP violation in the interference is present if the complex 
quantity A/, defined as 


Af = 


qAf 

q 

Af 

1L 

ii 

i 

-a 

Q 

p 

A f 


( 1 . 12 ) 


has a non-zero imaginary part. The right-hand expression is valid for a CP eigenstate / with eigen¬ 
value v\cp and 0 is the CP violating relative phase between q/p and Af/Af. The combination of 
CP violation in interference and in mixing is also called indirect CP violation. In general, it can 
be stated that CP symmetry is violated if A/ / 1. An excellent discussion on the different types of 
CP violation can be found in section 7.2.1 of Ref. [42]. As opposed to the strange and the beauty 
system, CP violation has not yet been discovered in the charm system. 

The occurrence of CP violation generally requires the interference of several amplitudes con¬ 
tributing to a particular process, which follows from writing a decay amplitude as 


A f = c(\+r f e l ^f + ^y 


(1.13) 


where the amplitude has been split in a leading and sub-leading part with a relative magnitude ry 
and relative strong and weak phases 8j and (j)f, respectively. Correspondingly, the CP-conjugate 
decay is given by 

Aj = c(\+r f e l{8 f -M), (1.14) 

where the weak phase changes sign and strong phase does not. Using the time-integrated decay 
rate 

= rr(M(t)->f)dt°c \Af\ 2 , 

Jo 


(1.15) 
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Figure 2: (Left) Tree-level diagram and (right) penguin diagram for a Z) () —> K decay. 


one obtains for the CP asymmetry 


r(M->/)+r(M->/) 


2 ry sin5/sin0y. 


(1.16) 


This equation shows that CP violation can only occur in the case that two (or more) amplitudes 
contribute to the decay (ry ^ 0) and that they have different strong and weak phases. 

The two amplitudes required for CP violation can for example be given by a tree-level decay 
and a penguin contribution (see Fig. 2). Penguin diagrams contain a qq pair of the same flavour, 
which reduces the number of decays to which they can contribute. For charm decays, these are the 
singly Cabibbo-suppressed c —> udd and c —> uss transitions, which occur for example in the decays 
D°—> n~n+ and D°^K-K + (see Sec. 3.2.2). 

Tree-level and penguin amplitudes generally have different strong phases and the different 
CKM elements present in the process lead to different weak phases, thus enabling CP violation. An 
alternative to two amplitudes interfering for the same decay exists for neutral mesons. In this case 
a final state can potentially be reached either through a direct decay or through the transition of 
the meson into its anti-meson followed by the anti-nreson decay to the same final state. Here, the 
difference in the weak phase can originate in the mixing process, while the meson and anti-nreson 
decay amplitudes may have different strong phases. An example of this scenario is the doubly 
Cabibbo-suppressed decay D° —> K + n~ (see Sec. 2.1.2 and Sec. 3.1.2). 


2. Charm mixing 

The studies of charm mesons have gained in momentum with the measurements of first evidence 
for meson anti-nreson mixing in neutral charm mesons in 2007 [21, 22]. Mixing is a process that 
changes the flavour quantum number by two units, e.g. charmness from +1 to — 1. The only flavour¬ 
changing interaction is the weak interaction, where the boson permits A F = 1 transitions. Thus, 
the only diagram for a A F = 2 process is a box diagram where two bosons are exchanged (see 
Fig. 3). This process is often referred to as short-distance effects in literature due to the small 
distance scales necessary for boson exchange. 

Mixing of D° mesons is the only mixing process where down-type quarks contribute to the box 
diagram. Unlike B-meson mixing, where the top-quark contribution dominates, the third generation 
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quark is of similar mass to the other down-type quarks. This leads to a combination of GIM 
cancellation [6] and CKM suppression [45, 46], which results in a strongly suppressed mixing 
process [47, 48, 49, 50, 51]. 

Another way for neutral meson mixing to proceed is through re-scattering effects. This can be 
seen as the decay of a neutral meson into a final state that is common to both flavours, followed 
by a recombination of the final-state particles to the anti-meson (see Fig. 3). These so-called long 
distance effects are likely to dominate charm meson mixing as they do not suffer from the same 
cancellations as the short distance effects. 

There are two approaches for theoretical calculations of charm mixing. The “inclusive” ap¬ 
proach is based an operator product expansion (OPE) in the ratio of the typical hadronic scale over 
the charm quark mass, A /m c [51, 52, 53, 54, 55, 56, 57, 58]. Due to the cancellations mentioned 
above it is higher order operators that give the largest contributions to the mixing parameters. Fur¬ 
thermore, it is not yet clear whether the expansion series really converges since A /m c ~ 0.25. Cal¬ 
culations of the charm meson lifetimes are being performed to test whether the OPE approach can 
properly reproduce the large difference between the D° and the D + lifetimes [58]. In the system, 
the OPE approach successfully predicted the width splitting of the two B® mass eigenstates [59] 
which has recently been confirmed by an LHCb measurement [60]. 

At this point it is useful to insert a comment on D lifetimes. The large ratio of the D + to the D° 
lifetime of about 2.5 indicates that the spectator quark plays an important role in the calculation of 
hadronic decay rates of charm mesons. In the OPE, terms depending on the spectator quark enter 
only at order (A/m c ) 3 underlining the importance of higher order calculations. The spectator quark 
is the only difference in the Cabibbo-favoured c —> sdu transitions. For D + decays the spectator 
quark, d, is the same as one occurring in the charm quark decay. Thus, there are two possible tree- 
level amplitudes, a colour-allowed and a colour-suppressed one, that can interfere. Their destructive 
interference leads to a large suppression of the partial widths of hadronic D + decays compared to 
their D° counterparts. 

Similarly to the decay width of a particle being given by the sum of all partial decay widths, 
the width splitting of the mass eigenstates can be calculated by considering the appropriate decay 
processes. The “exclusive” approach sums over intermediate hadronic states, taking input from 
models or experimental data [61, 62, 63, 64, 65, 66, 67, 68]. Also in this approach, different 
modes of the same St/(3) multiplet, e.g. D° —» hh' with h = K,n, lead to cancellations which is 
why their individual contributions have to be known to high precision. Due to the considerable 
mass of the D° meson, many different modes need to be taken into account simultaneously. Of 
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these, only phase space differences can be evaluated at the moment. Estimates indicate that mixing 
in the experimentally observed range is conceivable when taking into account SU (3)-breaking 
effects. However, neither the inclusive nor the exclusive approach have thus far permitted a precise 
theoretical calculation of charm mixing. A recent review of these calculations of charm mixing is 
given in Ref. [69]. 

It was discussed whether the measured size of the mixing parameters could be interpreted as 
a hint for physics beyond the standard model [70, 71, 72, 73, 74, 75, 76]. The biggest problem 
in answering this question is the non-existence of a precise standard model calculation. Effects of 
physics beyond the standard model were also searched for in numerous CP violation measurements 
and searches for rare decays both of which are covered in the remainder of this review. 

Given the smallness of the mixing parameters x and y in the charm system, it is instructive to 
expand the terms in Eq. 1.10 to obtain 


P(M°(t) —» M°) = P(M°(t) —>• M°) 
P(M°(t) —> M°) 

P(M°(t) —> M°) 


l e ~rt 

2 e 


2 - 


C 2 -y 2 ) 


(Ef) : 


„-r> [P+y 1 ) 
2 

,-Tt C 2 +y 2 ) 


(r;) 2 , 

(r t)\ 


(2.1) 


where terms of order (MT) 4 and (yrt) 4 have been neglected. Assuming a mixing parameter of the 
size of 0.5%, inspired by the latest world average [77], and a decay time corresponding to 20 D° 
lifetimes, this approximation corresponds to a relative error of 1%. 


2.1 Experimental measurements 

2.1.1 Flavour tagging 

Mixing of D° mesons can be measured in several different modes. Most require identifying the 
flavour of the D° at production as well as at the time of the decay. Tagging the flavour at production 
usually exploits the strong decay D* + —> D°k + (and charge conjugate 3 ) where the charge of the 
pion determines the flavour of the D°. The small amount of free energy in this decay leads to the 
difference in the reconstructed invariant mass of the D* + and the D°, dm = m D o K + —m D o, exhibiting 
a sharply peaking structure over a threshold function as background. While the small amount of 
energy transfer leads to a sharply peaking signal and thus to a high signal purity, the pion from 
the D* + decay has little kinetic energy and is thus commonly referred to as the soft pion. This is 
particularly important as low momentum particles can be swept outside detector acceptances by 
magnetic fields, thus reducing the overall detection efficiency. In addition, soft pions are more 
likely to be affected by detection asymmetries, which in the case of a tagging particle, can translate 
into a fake physics asymmetry. 

An alternative to using the D* + decay mode is tagging the D° flavour by reconstructing a 
flavour-specific decay of a B meson. At the LHCb experiment this approach is of interest due to 
differences in trigger efficiencies compensating for lower production rates. In particular the decay 
chain B-? D i] ji A, with X representing a non-reconstructed fraction of the final state, is very useful 
due to a high branching fraction of semi-leptonic decays and due to a very high trigger efficiency 

3 Charge conjugate decays are implicitly included henceforth. 
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for muons compared to hadrons. Another advantage of this tagging approach is that the efficiency 
for reconstructing the D° is independent of its flight distance with respect to its origin, while D° 
mesons originating directly in the pp interaction tend to have very low efficiency at low decay 
times. Thus, the two tagging methods not only provide complementary datasets but also improve 
the decay-time coverage, which improves the lever arm for measuring time-dependent effects. 

A third option available particularly at e + e“ colliders is the reconstruction of the opposite side 
charm meson in a flavour specific decay. This is based on the fact that quarks have to be produced in 
quark anti-quark pairs and thus every charm meson has to be accompanied by a hadron containing 
an anti-charm quark. For c e colliders operating at the threshold for producing y/(3770) states 
the correlation goes further and the final state are quantum-entangled D°-D° or D -D pairs. 

2.1.2 Mixing measurements 

Measuring the mixing rate: Theoretically, the most straight-forward mixing measurement is that 
of the rate of the forbidden decay D° —> K p v u which is only accessible through D°-D° mixing. 
The ratio of the time-integrated rate of these forbidden decays to their allowed counterparts, D°—» 
K p 1 v^, thus represents the rate at which D° mixing occurs. It is given, following Eq. 2.1, by 


/ 0 “T(Z>°(0-> K+ii-v^dt JoP(D°(t) - 

+ D°)\A f \ 2 dt 

R 

2 o | 2 

x~+y 

R 

fo°° r ( D °(t) K-p+v,)dt ./o°° P(D°(t) - 

A D°)\A f \ 2 dt " 

P 

2 

P 


where |A/| 2 = \Aj\ 2 are the decay rates of the decay D° —> K p and its CP-conjugate, re¬ 
spectively, which cancel in the ratio. Furthermore, the mixing rate R m = (x 2 ±y 2 )/2 has been 
introduced. 

As this requires very large samples of D° mesons no measurement has thus far reached suffi¬ 
cient sensitivity to see evidence for D° mixing. The most sensitive measurement to date has been 
made by the Belle collaboration [78] to P m = (1.3 ±2.2 ±2.0) x 10 4 , where the first uncertainty is 
of statistical and the second is of systematic nature 4 . Using the current world average values [77], 
this can be compared to the expected rate of approximately R m = (0.3 ± 0.1) x 10 4 . 


Mixing through interference: Related to the semi-leptonic decay is the decay D° —> K ’ ji . 
This decay can proceed either via mixing followed by the Cabibbo-favoured (CF) decay D° —> 
K + n~ or via the direct decay D° —> Kn\ which is doubly Cabibbo-suppressed (DCS). These 
two amplitudes interfere and the combined process is referred to as wrong-sign (WS) decay (see 
Fig. 4). The corresponding Cabibbo-favoured process is calles right-sign (RS) decay. This leads to 
the WS decay rate, using Eq. 1.8, 


F (D°(t)^K + n~) 


~f- M/ ± /+ (?)Ay 


2 


(2.3) 


where A j is the DCS amplitude of the D° decay and A j is the CF amplitude of the D° decay. 
Assuming CP conservation, the time-dependent decay rate of the WS decay is 


r(D°(f) -A K + TI~) = e~ Tt (R d ± VR^y'rt + R m (rt) 2 ) , (2.4) 

4 This notation is applied to all results where two uncertainties are quoted. 
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Figure 4: (Left) Schematic of the amplitudes contributing to RS and WS decays. The bold lines 

indicate favoured amplitudes, while the thin connections represent suppressed amplitudes. (Right) Decay¬ 
time distribution of the ratio of WS to RS D° —> K 1 7T f decays leading to the first single-measurement 
observation of charm mixing. Reproduced from Ref. [80]. 


where Rn = \Aj/Aj\ 2 is the ratio of the DCS to the CF rate. The mixing parameters are rotated by 
the strong phase difference between the DCS and the CF amplitude, 8kk, leading to the observable 
y' = y cos 8 Kn - jc si n 8 Kn [79]. 

First evidence for charm mixing was obtained by the BaBar collaboration in March 2007 
based on a measurement of wrong-sign D° —> K + n~ decays [21]. While the two observables x' 2 
and y' were not measured to be non-zero with a significance exceeding three standard deviations, 
the no-mixing point x' 2 = / = 0 is excluded with a significance of 3.9 standard deviations. Fur¬ 
ther supporting evidence was also obtained in the same measurement by the CDF collaboration 
later in the same year [81]. It then took until late 2012 for the LHCb collaboration to perform a 
measurement excluding the no-mixing hypothesis by 9.1 standard deviations, thus providing first 
single-measurement observation of mixing in the charm system [80]. This result can be seen in 
Fig. 4 where the clearly non-zero slope indicates the presence of mixing. Subsequent updates led 
to the following set of latest measurements 


x' 2 in 10 3 / in 10 3 

CDF [82] 0.08±0.18 4.3±4.3 

LHCb [83] 0.06 ±0.05 4.8 ±1.0 

Belle [84] 0.09 ±0.22 4.6 ±3.4 

Similarly, the CF and DCS amplitudes can also lead to excited states of the same quark content. The 
decay D° — > K + n~n° is the final state of several such resonances, e.g. K + p~ and K* + n~ . Thus, 
by studying the decay-time dependence of the various resonances a mixing measurement can be 
obtained. The BaBar collaboration achieved a measurement showing evidence for D° mixing [85] 
with central values of x" = (26.1+^ ± 3.9) x 10~ 3 and y" = (-0.6+^ ± 3.4) x 10~ 3 , where the 
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rotation between the observables and the system of mixing parameters is given by a strong phase 
difference as 


a xcos §K+ft- jjO Ty sin 8 ^+^- (2-5) 
y" =y cos 5 K+n - n o — xsin5 K+K - K a. (2.6) 

The significant advantage of this analysis over that using two-body final states is that both mixing 
parameters are measured at first order rather than one at first and one at second order. 

The strong phase differences are not accessible in these measurements but have to come from 
measurements performed using quantum-correlated D-D° pairs produced at threshold. Such mea¬ 
surements are available from CLEO [86, 87, 88, 89] and BESIII [90]. In addition, further con¬ 
straints on these strong phase differences can be obtained from the combination of several mea¬ 
surements that share the underlying mixing parameters but are subject to different strong phase 
differences. This can be for example the combination of D° —> K '• 71 and D° —> K l 71 n° together 
with other measurements that are not affected by strong phases (see below). 


Mixing in CP eigenstates: In the absence of CP violation the physical eigenstates are CP eigen¬ 
states and therefore the width difference can be accessed directly through a measurement of the 
effective lifetime of a CP eigenstate with respect to the lifetime of a flavour-specific state, which is 
not affected by mixing, using the observable 


yep = 


Tcpa 

r 


r 2 .i Ar 

~ 1 “ T” ~ 1 = 2 f = y ‘ 


(2.7) 


In a more general scenario allowing also for CP violation the effective lifetimes of D° (D°) mesons 
into final states that are CP eigenstates, f (f), lead to 


yep 


2 t d 

T + f 


-1 


flCP 



ycos (j) 


. , 

—isinffl , 
2 r 


( 2 . 8 ) 


where To is the D° lifetime and 


t Mm — 




(2.9) 


which leads to \q/p\ ±2 ~ I ± a m [91]. As the CP-violating contributions a m and 0 enter only 
at second order, measurements of yep are among the most powerful constraints of the mixing 
parameter y. 

Following only a few weeks after the BaBar measurement using D° —> K n decays, a mea¬ 
surement of the Belle collaboration of the parameter yep was among the first to provide evidence 
for mixing in the charm system [22]. The latest results are 
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_ yep in 10~ 3 

LHCb [92] 5.5 ±6.3 ±4.1 

BaBar [93] 7.2 ±1.8 ±1.2 

Belle [94] 11.1 ±2.2±1.1 

The BaBar collaboration has added the larger sample of untagged events to their analysis; 
however, with limited gain in sensitivity due to larger systematic uncertainties for the untagged 
sample which has lower purity compared to the D* -tagged events. The LHCb result is based on 
only about 1 % of their total dataset and significant improvements in precision are expected from 
future updates. 

The central values of these results yield an average value of yep, which is greater than that of 
the mixing parameter y. CP violation can only mildly increase yep over y (see Eq. 2.8) and thus 
such a tension would not easily be explained. However, based on the latest results, this discrep¬ 
ancy is at the level of two standard deviations even for the Belle result and thus considered not 
significant [95]. 

Another possibility of measuring yep is using the decay mode D° —> K^K K 1 . The Belle 
collaboration have published a measurement in which they compare the effective lifetime around 
the (j) resonance with that measured in sidebands of the K K + invariant mass [96]. The effective 
CP eigenstate content in these regions is determined with two different models. Their result is 
y C p = (1.1 ±6.1 ±5.2) x 1(T 3 . 

Similarly, one can measure yep in a phase-space integrated way provided that the decay is 
close to being a CP eigenstate [97]. Based on measurements of the CLEO collaboration the decay 
D° —> K n n {) has been confirmed to be practically CP even [98]. Other potential candidates 
are D° —> n n n JT and D° —y n 1 n { \ for which the effective CP content remains to be 
measured. 

Direct measurements of the mixing observables: The decay D° —> K^K~K + and more so the 
decay D° — > K$ 7l~7t + give excellent access to the mixing parameters x and y individually. This is 
achieved through the simultaneous measurement of the decay-time evolution and resonance ampli¬ 
tudes in the Dalitz plot. At the same time measurements of these final states allow a determination 
of parameters of indirect CP violation as discussed in the following section. Under the assumption 
of no CP violation Belle and BaBar have measured 

x in 10~ 3 y in 10~ 3 

BaBar [99] 1.6±2.3± 1.2±0.8 5.7±2.0± 1.3±0.7 

Belle [100] 5.6± 1.9+°; 3 +°;® 3.0± 1.5+°j+°; 3 

where the last uncertainty in each measurement is a model uncertainty. With this analysis being 
arguably the most complicated charm analysis, there is no result from the LHCb collaboration to 
date. However, their large dataset should enable significant improvements of this measurement, 
which will be of particular benefit to contraining the mixing parameter x. 
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Combined results: As discussed above, no single measurement has thus far been able to pin down 
both mixing parameters x and y. Therefore, all available measurements are being combined by the 
Heavy Flavor Averaging group in order to obtain maximal sensitivity to the underlying parameters. 
In doing so, measurements of strong phase differences are used, but also the combination itself is 
overconstraining and thus helps to pin down these phases. 

Under the assumption of no CP violation in mixing or decays, the world average of the mixing 
parameters is x = (4.9^} 5 ) x 1(U 3 and y = (6.2 ±0.8) x 1(U 3 [77]. This result affirms the conclu¬ 
sion that the CP even final state is the shorter-lived. However, the periodic oscillation of D° and D° 
mesons and vice versa, which is governed by the parameter x, has not yet been established. The 
no-mixing hypothesis, i.e. R m = 0, is excluded by well over 12 standard deviations. 

2.1.3 Measurement techniques 

Although all measurements discussed above are based on the decay-time evolution of D° decays, 
not all of the have been performed by a fit to an exponential distribution or modulation thereof. 
Essentiall all measurements rely on the comparison of the decay-time evolution of two related 
processes. The measurement based on semi-leptonic decays can be performed time-integrated as 
its time dependence is given by a single term proportional to t 2 . 

The measurement based on wrong-sign D° —y K ' n decays and that of yep is based on the 
comparison of the time evolution of two related decay modes. This can be done by fitting the 
appropriate model, an exponential multiplied by a second order polynomial for the WS K + n 
decay or a simple exponential for CP eigenstates, and comparing that to the lifetime measured in 
the CF D° —> K 71 1 decay. An alternative is to divide the data into bins of decay time and to 
measure the ratio of signal yields in these bins. This ratio can then be fitted using a second order 
polynomial for the WS K + n decay or a linear function to extract yep. 

Both techniques have been used in the different measurements and sometimes they were taken 
to provide mutual cross-checks in the same measurement. A major advantage of the binned ratio 
method is that it cancels decay-time acceptance effects when they can be assumed to be sufficiently 
similar for the two modes being compared. This is generally the case for the WS K + n decay 
measurement, but not for yep as the latter compared decays with different final state particles. The 
decay-time acceptance argument is particularly relevant to hadron colliders where signal candidates 
are usually required to originate from some distance away from the primary collision point. At e 
e + colliders the D° mesons are generally less boosted and therefore the decay-time resolution is 
worse. In this context bin migration becomes an important factor to take into account, while the 
resolution is readily taken into account in an exponential fit. 

The situation is more complicated for multi-body final states. These generally require an 
analysis of the decay-time dependence of the different contributions to the phase space. This can 
be obtained by two different groups of methods. The first uses a model of the phase space and fits 
a dedicated model to the decay-time distribution each of the resonances or groups of resonances 
contributing to the phase-space model. The decay-time models incorporate the common underlying 
theory parameters depending on the amplitudes contributing to the specific resonance. 

The alternative method is based on measuring the decay-time dependence in specific regions 
of phase space. These regions were chosen to enclose areas of similar strong phase difference and 
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the choice of these regions is inspired by a phase-space model. The effective CP content of each 
of these regions needs to be known, which can be obtained through measurements of quantum- 
entangled D° mesons, such as in Ref. [101]. Using this information, the decay-time distributions 
can be linked to the underlying theory parameters once more. The challenge of the former approach 
lies in assigning an uncertainty to the accuracy of the phase-space model, while the latter depends 
on accurate measurements of the CP content of these regions. While the accuracy of models may 
improve in the future through a combined effort of experiment and theory, any improvement on the 
second approach currently relies on further measurements made by the BESIII collaboration. 


3. Charm CP violation 


3.1 Indirect CP violation 


3.1.1 Measuring indirect CP violation directly 

Indirect CP violation is the deviation from unity of q/p, which comes with the fact that the mass 
eigenstates |Di) and /U) are no longer CP eigenstates. This is measured most straight-forwardly 
by comparing the decay-time dependence of D° and D° decays to CP eigenstates. Writing the 
effective lifetimes of the D° and D° decays as f and t, respectively, this leads to the observable 


A r = I —\ ~ rj C p 

Z+Z 


- (a m + fl , d)vcos 0 — *sin 0 


_ + 2 a< *y cp > 


(3.1) 


which has contributions from both direct and indirect CP violation [91, 102]. 

All available measurements of Ap to date have been based on the two-body CP eigenstates 
K K and n X r . The BaBar and Belle collaborations have published measurements based on 
the combination of these final states. They obtain Ap = (0.9 ± 2.6 ± 0.6) x 10 ~ 3 [93] and Ap = 
(—0.3 ±2.0 ±0.8) x 10 3 [94], respectively. The LHCb collaboration has measured Ap separately 
for the two decay modes and separately for two flavour-tagging approaches. The results based 
on a sample where D* + —y D°n + decays are used for tagging the D° flavour are A\ (K K 1 ) = 
(—0.4±0.6±0.1) x 10 ~ 3 and Ay{K~K + ) = (0.3± 1.1 ±0.1) x 10 ~ 3 [103], representing the most 
precise measurement of a charm CP asymmetry to date. As this result is based on only one third 
of LHCb’s integrated luminosity, further improvement may be expected with existing data. In 
addition, a recent analysis is based on a complementary data sample where the D° flavour is tagged 
using the charge of the muon from semi-leptonic b-hadron decays involving D° mesons. The results 
areA r (K-K + ) = (-1.3±0.8±0.3) x 10 ^ 3 andA r (7W£ + ) = (-0.9± 1.5±0.3) x 10 ' 3 [104], All 
individual measurements are in good agreement with CP conservation. A naive average assuming 
negligible final state differences yields Ap = (—0.6 ± 0.4) x 10~ 3 . 

Using current experimental bounds values of Ap up to C/{ 10 4 ) are expected from theory [102, 
105]. It has however been shown that enhancements up to about one order of magnitude are pos¬ 
sible, for example in the presence of a fourth generation of quarks [51] or in a little Higgs model 
with T-parity [105]. The current level of precision starts to put bounds on this parameter space and 
can be expected to be further improved with the analysis of the full LHC run-1 dataset of the LHCb 
experiment. 

Eventually, the interpretation of the observable Ap requires precise knowledge of both mixing 
and CP violation parameters. The relative sensitivity to the CP-violating quantities in Ap is limited 
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by the relative uncertainty of the mixing parameters. Therefore, to establish the nature of a potential 
non-zero measurement of Ap it is mandatory to have measured the mixing parameters with a relative 
precision of about 10%. 


3.1.2 CP asymmetries in D° —tK % decays 

Further powerful access to CP asymmetries lies in D° —> K + n~ decays. When conducting the 
measurement described in Sec. 2.1.2 separately for D° and D° decays one can define the set of 
observables Rx a± , and y !± , where the ± refers to the parameter for the D° and D° decay, respec¬ 
tively. Following from Eq. 2.3, these are related to the underlying parameters by 




/± 


y 


i± 


q 

p 

q 

p 


±t 

{x' cos (p ±/sin0) 
±t 

(/ cos (p =p x sin (p ), 


(3.2) 


and a non-zero difference between R ^ and R [} would indicate CP violation in the DCS decay. For 
small phases (p this gives direct access to the magnitude of | q/p\ and thus CP violation in mixing. 
For larger weak phases the sensitivity is worse since the sensitivity to x' is much worse compared 
to /. Several measurements of these parameters have been made by Belle [106], BaBar [21], and 
LFICb [83] and none of them shows any indication of CP violation. 


3.1.3 Direct measurements of | q/p\ and <f> 

The analyses of the decays D° —> K^n n' and D °—> K^K K l offer separate access to the pa¬ 
rameters x, y, \q/p\ and arg (q/p) and are one of the most promising ways of obtaining precise 
mixing measurements. These analyses require the determination of the decay-time dependence of 
the phase space structure as described for the mixing measurements in Sec. 2.1.2. The only existing 
measurement to include the CP violating parameters has been performed by the Belle collabora¬ 
tion to yield \q/p\ = 0.90+q 04^0 05 anc * 0 = ( — 6± 11 ±3l4) 0 [100]. Other measurements 
were performed by the CLEO [107] and BaBar [99] collaborations assuming CP conservation and 
thus extracting only x and y. With the data samples available and being recorded at LFICb and 
those expected at future flavour factories, these measurements will be very important to under¬ 
stand charm mixing and CP violation. However, in order to avoid systematic limitations it will be 
important to reduce model uncertainties or to improve model-independent strong-phase difference 
measurements which are possible at BESIII. 


3.1.4 Combining indirect CP violation measurements 

Direct access to CP violation parameters is limited to the one class of measurements described in 
the previous paragraph. Therefore, extracting the underlying theory parameters from a combined 
fit to all relevant results is required. Figure 5 shows how different measurements contribute to the 
current world average as computed by the Heavy Flavor Averaging Group. The brown contour in 
the background essentially combines measurements of D° —> K { ^n which lead to an elliptical 
constraint, and measurements of Ap, which constrain a diagonal in this parameter space. The ma¬ 
genta contour shows the power of measurements based on D° —> K ' n decays alone, underlining 
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Figure 5: Contributions of the combination of different measurements to the extraction of the parameters 
\q/p\ and tj). Reproduced from Ref. [95]. 


the particular sensitivity to \q/p\ for small values of (j). Finally, the green contour in the foreground 
is the full combination of all measurements. The official world average is \q/p\ = 0.93 1 and 
0 = (-8.7t|;l)° [77], 

Under the assumption that there is no CP violation in the corresponding decay rates one can 
construct the following relationship [71, 102, 108]: 


tan (j) 


i-lff/pl 2 x 

i + \q/p \ 2 y 


(3.3) 


This equation can be used to reduce the number of mixing and indirect CP violation parameters 
from four to three. Using this to remove 0 the world average on \q/p\ becomes 1.007 and 
correspondingly, removing \q/p\ one obtains (j) = (—0.30+q 5q)° [77], which shows an impressive 
improvement in sensitivity. 


3.2 Direct CP violation 


Direct CP violation is the asymmetry of a decay rate to its CP conjugate, defined as 


aa>(D-tf) 


r(D->/)-r(D->/) 

r(D->f)+r(D->f)’ 


(3.4) 


for the decay of a D meson to a final state /. Direct CP asymmetries are generally measured in 
decay-time integrated measurements. For neutral mesons, the decay-time distribution of the data 
has to be taken into account to estimate the contribution from indirect CP violation as will be further 
explained in Sec. 3.3. 


3.2.1 Experimental challenges 

Experimentally, there are a number of effects that can generate asymmetries and thus obscure 
access to the physics asymmetry acp■ Under the assumption of small asymmetries, these effects 
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are additive to first order and can be written as 

a raw = a CP ± a det,f ± ^det.tag ± ^prod; (3-5) 

where fldet.f and fldet.tag are detection asymmetries of the final state and the tagging process, respec¬ 
tively, and where a proc | is the production asymmetry of the initial state. 

The production asymmetry is only an issue at matter anti-matter asymmetric machines. These 
are pp colliders like the LHC or fixed target experiments where a matter beam is shot on a matter 
target. However, pp colliders such as the Tevatron can have a forward-backward asymmetry as 
protons and anti-protons have a preferred direction. A forward-backward asymmetry exists also at 
e + e~ colliders where its origin lies in the interference of particle production via virtual photons and 
via virtual Z bosons. Experimentally, production asymmetries are measured using control modes 
with the same initial state as the signal decay. 

The detection asymmetries have two origins. These are detector asymmetries where the detec¬ 
tor geometry, mostly in combination with a magnetic field, leads to a higher probability of detecting 
particles from one charge compared to the other. One example can be caused by staggering of de¬ 
tector elements, which can lead to an acceptance gap for particles travelling in a specific direction 
while it is fine for particles traversing the detector at a different angle. The second component are 
interaction asymmetries, which are based on the fact that particles of one charge have a different 
probability of interacting with the detector than the particle of the other charge. One example are 
K mesons for which the strange quark can produce hyperons and the anti-up quark can annihilate 
with the detector material, whereas the anti-strange quark of the K + can neither produce hyperons 
or annihilate with quarks in the detector material. Detection asymmetries may intrinsically cancel 
if two particles of the same kind and with opposite charge are present in the final state, such as in 
the decay D° —y K K + . Otherwise, detection asymmetries need to be measured in control modes 
and subtracted off. 

Due to the underlying mechanisms that cause production and detection asymmetries they may 
vary with the particle’s kinematics or its position in the detector. Therefore, any subtraction of 
asymmetries needs to ensure sufficient overlap in the relevant quantities to guarantee proper can¬ 
cellation. 

3.2.2 Experimental measurements 

CP violation in two-body final states: Keeping in mind the above, a very powerful observable is 
the difference in CP asymmetries of the decays of D° mesons into two charged pions or kaons. In 
this case all nuisance asymmetries cancel to first order and the measured difference is equal to the 
difference in physics observables, defined as 


— acp(K K ) dcp(7l 71 ) — 6! raw (K K ) ^raw(^t 71 ). 

An initial measurement by the LHCb collaboration was the origin of an increase in experi¬ 
mental and particularly theoretical activity in this field. The result of A cicp = (—8.1 ±2.1 ± 1.1) x 
10 3 [109] provided first evidence for CP violation in the charm sector. Subsequent measurements 
did thus far not confirm this evidence and the latest world average of (—2.5 ± 1.0) x 10 3 already 
excludes percent level CP violation. Measurements from BaBar [110], Belle [111], CDF [112], and 
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LHCb [113] are based on D* ' -tagged decays, while the LHCb collaboration has also published a 
measurement based on D° mesons originating in semi-leptonic b-hadron decays [114]. 

The observable A acp gives access to the difference in direct CP violation of the two decay 
modes through 



(3.6) 


where x is the nominal D° lifetime, X = (X(K~K + ) + X{n^ K + ))/2, and AX = X(K~K + ) — 
X(n^iz + ) for X = (a^p, (t) ,Ap) [91]. Here, Ar is multiplied with the difference of mean decay 
times as indirect CP violation leads to different effective lifetimes. Thus, its impact on a time- 
integrated asymmetry depends on the decay-time distribution of the measured data sample. Equa¬ 
tion (3.6) assumes the CP-violating phase 0 to be universal. For a small non-zero difference in this 
phase between the two final states, A (f>f, an additional term of the form xA<j>f(t) / X arises as pointed 
out in Ref. [102]. Given a typical variation of (t) / x between 1 and 2.5 for the different experiments 
the contribution of A (j)f is suppressed by x(t)/x ~ 10 ~ 2 . 

An additional beneficial feature of the observable A acp is that the two contributing asymme¬ 
tries are expected to be roughly of equal magnitude and opposite sign, due to the CKM structure of 
the decays. Thus, A acp would measure twice the magnitude of one of the individual asymmetries. 
Nevertheless, in addition to measuring the difference, the knowledge of the individual asymmetries 
is required to pin down the source of potential CP violation. Several measurements of the individual 
asymmetries exist, with the most recent and most precise one by the LHCb collaboration [114]. 

The challenge in measuring these individual asymmetries clearly lies in controlling all nui¬ 
sance asymmetries to extract the acp from the measured asymmetry. The LHCb measurement is 
based on flavour-tagging through semi-leptonic b-hadron decays and their production asymmetry 
is controlled through a measurement of Cabibbo-favoured D° —> K n 1 decays produced in the 
same way. While this also cancels the detection asymmetry of the tagging muon, it introduces a 
detection asymmetry due to the K n + final state. This in turn is measured in the difference of the 
CF decays D + —> K n n~ and D + —> K^n , where in the latter decay the Kg decays into 7i + n~. 
This leaves the asymmetry introduced by the Kg due to regeneration in the detector material and 
due to CP violation in kaons. Both effects are well known and can be calculated and subtracted to 
sufficient precision. The latest world average is compatible with both asymmetries being of equal 
magnitude and opposite sign with a Kn = (1.4 ± 1.5) x 10 3 and a nn = (—1.1+J 3 ) x 10 ~ 3 [77]. 

CP violation in two-body final states of charged D mesons: Decays of D + and D+ into a Kg and 
either a K + or a n + are closely related to their D° counterparts. Measurements of time-integrated 
asymmetries in these decays are expected to exhibit a contribution from CP violation in the kaon 
system. As has been pointed out [115] this contribution depends on the decay-time acceptance of 
the Kg. This can lead to different expected values for different experiments. For Belle [116], the 
expected level of asymmetry due to CP violation in K° —K° mixing is —3.5 x 1 0 3 . For LHCb on 
the other hand, there is no significant asymmetry induced by kaon CP violation [117] as the LHCb 
acceptance, for Kg reconstructed in the vertex detector, corresponds to about 10 % of a Kg lifetime. 
CP violation searches in the decays D + —>• 7^71 + [118, 119, 120] and D+ —y Kg7 r + [1 19, 120] show 
significant asymmetries. However, these asymmetries are fully accommodated in the expected CP 
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violation of the kaon system. These measurements do not show any hint for an asymmetry in D 
decay amplitudes. This is confirmed by recent LHCb measurements in which kaon CP violation is 
negligible [121]. 


CP violation in multi-body final states: Another group of channels suitable for CP violation 
searches is that of decays of D + and D+ mesons into three charged hadrons, namely pions or kaons. 
Here, CP violation can occur in quasi two-body resonances contributing to these decay amplitudes. 
Asymmetries in the Dalitz-plot substructure can be measured using an amplitude model or using 
model-independent statistical analyses [122, 123, 124, 125, 126]. The latter methods allow CP 
violation to be discovered while eventually a model-dependent analysis is required to identify its 
source. The two types of model-independent analyses differ in being either binned [122, 123, 124] 
or unbinned [124, 125, 126] in the Dalitz plane. 



Figure 6: (Left) Dalitz plot of D + —> n K K ' candidates and (right) asymmetry significance Sep for the 
same candidates. Reproduced from Ref. [124]. 


The binned approach computes a local per-bin asymmetry significance, Sep, and judges the 
presence of CP violation by the compatibility of the distribution of local asymmetries across the 
Dalitz plane with a normal distribution. An example of a Dalitz plot and corresponding Scp distri¬ 
bution is shown in Fig. 6. This method obviously relies on the optimal choice of bins. Bins ranging 
across resonances can lead to the cancellation of real asymmetries within a bin. Too fine binning 
can reach the limit of statistical sensitivity, whereas too coarse binning can wash out CP violation 
effects by combining regions of opposite asymmetry. A model-inspired choice of binning is clearly 
useful and this does not create a model-dependence in the way that fitting resonances directly does. 
This method does not yield an easy-to-interpret quantitative result. This issue has been discussed 
in a recent update of the procedure [127]. 

The unbinned asymmetry search calculates a test statistic that allows the assignment of a p- 
value when comparing to the distribution of the statistic for many random permutations of the 
events among the particle and anti-particle datasets [125]. Moreover, being unbinned, there is 
no need for a model-inspired choice of binning. The drawback of this method is its requirement 
on computing power. The calculation of the test statistic scales as the square of the number of 
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Figure 7: Local asymmetry significances in the Dalitz plot of Z)°—>■ n K 1 candidates. Reproduced from 
Ref. [126], 


events. In a first implementation of the method this challenge has been overcome by the usage 
of graphics processing units (GPUs) [126]. Figure 7 shows the distribution of local asymmetry 
significances contributing to the above mentioned p-value calculation. The unbinned technique 
exploited in Ref. [124] can be seen as a simplification of the above procedure since it restricts the 
comparison to events that are nearest neighbours in phase space. This comparison being applied to 
the order of 100 nearest neighbours for millions of signal candidates clearly reduces the computing 
requirement, but similarly reduces the method’s sensitivity. 

Beyond three-body final states similar analyses can be performed in decays into four hadrons, 
e.g. decays of D° into four charged hadrons. This too gives access to interesting resonance struc¬ 
tures that may exhibit significant CP asymmetries. However, rather than having a two-dimensional 
Dalitz plane, the phase space for four-body decays is five-dimensional (see e.g. Ref. [128]). This 
poses not only a challenge on the visualisation but also on any binned approach due to rapidly 
decreasing sample sizes per bin. Also, the phase-space substructure can no longer be described 
only by interfering amplitudes of pseudo two-body decays as also three-body decays may con¬ 
tribute. The LHCb collaboration has released a first model-independent search for CP violation 
in the decays D° —> K K' n' % and D° —y n n l Ji n without finding any hint of CP non¬ 
conservation [129]. 

Theory of charm CP violation: The initially large central value for A acp led to an intensive 
discussion on the potential origin of this CP asymmetry. In particular, it could not be ruled out 
that such a value could be caused by standard model effects. On the other hand, since the size of 
standard model effects could not be pinned down precisely it could similarly not be ruled out that 
at least part of the effect was caused by physics beyond the standard model. With the current world 
average being significantly smaller, it is more conceivable that an asymmetry of this size is in the 
range of standard model effects but other sources can still not be ruled out. 

Within the standard model the central value can only be explained by significantly enhanced 
penguin amplitudes. This enhancement is conceivable when estimating flavour SU{ 3) or U-spin 
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breaking effects from fits to data of D decays into two pseudo scalars [130,131, 132, 133, 134,135]. 
U-spin describes a subgroup of SU (3) that exchanges d-quarks and s-quarks. However, attempts of 
estimating the long distance penguin contractions directly have failed to yield conclusive results to 
explain the enhancement. 

Lattice QCD has the potential of assessing the penguin enhancement directly. However, sev¬ 
eral challenges arise which make these calculations impossible at the moment [136, 137, 138, 139, 
140, 141, 142]. Following promising results on K —> Tin decays, additional hurdles arise in the 
charm sector as nn and KK states mix with rj r \, 4n, 6n and other states. Possible methods have 
been proposed and first results may be expected within the next decade. 

General considerations on the possibility of interpreting A acp in models beyond the standard 
model have lead to the conclusion that an enhanced chromomagnetic dipole operator is required. 
These operators can be accommodated in minimal supersymmetric models with non-zero left-right 
up-type squark mixing contributions or, similarly, in warped extra dimensional models [143, 144, 
145, 146, 147, 148, 149]. Tests of these interpretations beyond the standard model are in the 
focus of ongoing searches. One promising group of channels are radiative charm decays where the 
link between the chromomagnetic and the electromagnetic dipole operator leads to predictions of 
enhanced CP asymmetries of several percent [150]. 

Another, complementary, test is to search for contributions beyond the standard model in A7 = 
3/2 amplitudes. This class of amplitudes leads to several isospin relations which can be tested in a 
range of decay modes, e.g. D -» Tin, D —> pn, D —> KK [130, 151]. Several of these measurements, 
such as the Dalitz plot analysis of the decay D° —>• n n 7 t (l , have been performed by BaBar and 
Belle and will be possible at LHCb as well as future e + e~ machines. 

The Belle collaboration has recently published measurements of the challenging decays to 
neutral final states D° —» 7i°7i° and D°— > Kg7Z°. They find A(7T°7 z: 0 ) = (—0.3 ± 6.4 ± 1.0) x 10 ~ 3 
and A(^71°) = (-2.1 ± 1.6 ±0.7) x 10 " 3 [152], 

Beyond charm physics, the chromomagnetic dipole operators would affect the neutron and 
nuclear electric dipole moments (EDMs), which are expected to be close to the current experimental 
bound [145]. Similarly, rare flavour-changing neutral current (FCNC) top decays are expected to be 
enhanced. Furthermore, quark compositness can be related to the A acp measurement and tested in 
dijet searches. Current results favour the new physics contribution to be located in the D° —> K K 1 
decay as the strange quark compositness scale is less well constrained [153]. 

In the light of the recent measurements it is evident that there are four directions to pursue: 

• More precise measurements of A acp and the corresponding individual asymmetries are re¬ 
quired to confirm the size of potential CP violation in these decays. 

• Further searches for time-integrated CP violation need to be carried out in a large range of 
decay modes that allow the identification of the source of the CP asymmetry. 

• Searches for time-dependent CP asymmetries, particularly via more precise measurements 
of Ar, D° —y K 1 ji decays, and D° —>• K^k n decays are required to establish or constrain 
the size of indirect CP violation. 

• And finally a more precise determination of the mixing parameters, in particular x, is required 
to fully exploit the other observables. 
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3.3 Interplay of mixing, direct and indirect CP violation 

Following Eqs. (3.1) and (3.6) it is obvious that both Ap and A acp share the underlying CP-violating 
parameters. Allowing for a non-universal CP-violating phase 0 one can write 

A r(f) = ~a cp - a cp(f)ycp ~ (3-7) 

a CP (f)=a%(f)-A r (f)^, (3.8) 

Ac if", Alp—- ,\|——, (3.9) 

T t 

where again X = (. X(K~K + ) +X(7i~7i + ))/2 and AX =X(K^K + ) —X(ji~n + ) forX = (afp. (t),Ar). 
It is expected that, at least within the standard model, one has afj,(K~K + ) = — a^,(n~7C + ) and thus 
A\- = —Cpp. This set of equations shows that it is essential to measure both time-dependent (Ap) 
and time-integrated asymmetries (acp) separately in the decay modes D° —> K K + and D° —y 71 71 1 
in order to distinguish the various possible sources of CP violation. Currently, the experimental 
precision on Ar is such that there is no sensitivity to differences in the contributions from di¬ 
rect CP violation to measurements using K K + or 71 n final states. Hence, the approximation 
Ap = Ap ~ Ar(K~K + ) « Ap( 7 z;~ 7 r + ) can be used to obtain 

A r = —cicp (3-10) 

Aa CP = Aafp fl+yci^l + af^. (3.11) 

These equations have been used by HFAG to prepare a fit of the direct and indirect CP violation 
contributions [77] as shown in Fig. 8 . This fit yields a confidence level of about 1.8 x 10 2 for 
the no CP violation hypothesis with best fit values of A a^p = (—2.57 ± 1.04) x 10 1 and a' ( "p = 
(0.58 ± 0.40) x 10~ 3 . While there may be direct CP violation at the 10 -3 level, indirect CP violation 
is constrained to be less than 10 -3 . 

The ultimate goal of mixing and CP violation measurements in the charm sector is to reach 
precisions at or below the standard model predictions. In some cases this requires measurements in 
several decay modes in order to distinguish enhanced contributions of higher order standard model 
diagrams from effects caused by new particles. 

A dominant contribution to indirect CP violation measurements is given by the observable Ar 
(see Eq. (3.7)). The CP violating parameters in this observable are multiplied by the mixing param¬ 
eters x and y, respectively. Hence, the relative precision on the CP violating parameters is limited 
by the relative precision of the mixing parameters. Therefore, aiming at a relative precision below 
10% for the CP violation quantities and taking into account the current mixing parameter world av¬ 
erages, the target precision for the mixing parameters is 2 — 3 x 10 4 , corresponding to about 5% of 
the current uncertainty of the world averages. With standard model indirect CP violation expected 
of the order of 10 4 , the direct CP violation parameter contributing to Ap has to be measured to an 
absolute precision of 10 3 in order to distinguish the two types of CP violation in Ap. 

Expectations for direct CP violation vary for different decay modes. In addition to the modes 
discussed so far, radiative decays have been pointed out as promising probes for physics beyond 
the standard model. Possible CP asymmetries exceeding 1% are expected in D° —> Vy (V = p,(j)) 
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Figure 8: Fit of A a^p and Reproduced from Ref. [77]. 

decays in the tails of the invariant mass distribution of the vector resonance [150]. These regions 
have particularly high sensitivity to contributions through dipole operators. 

For multibody final states the aim is clearly the understanding of CP asymmetries in the in¬ 
terfering resonances rather than global asymmetries. Of highest interest are those resonances that 
are closely related to the two-body modes used in A acp, for example the vector-pseudoscalar reso¬ 
nances K*K and pn. The measurement of further suppressed resonances is of interest as well since 
those have no contributions from gluonic penguin diagrams, thus allowing to constrain the source 
of CP violating effects. 

3.4 CP violation through triple product asymmetries 

In addition to rate asymmetries, whether global or local to a region of phase space, decays 
with a phase space, which is described by more than two degrees of freedom, can be used to study 
asymmetries in triple products. Examples for triple products are (see e.g. Ref. [2]) 

P\ ■ (pi x pi), s-[p\ x p 2 ), (3.12) 

where pi are the momenta of the decay products of e.g. a four-body decay and ?is a spin vector. 

Such triple products are odd under parity transformation as well as under time reversal. In 
the literature, these triple products are often referred to as T -odd variables, which indeed they are. 
Others prefer just to refer to them as P-odd variables, since only decays are analysed and not their 
time reversal. A detailed dissection of triple-products and derivative observables can be found in 
Ref. [154], 
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In order to search for CP violation, the asymmetry for the triple product Cp 

_ r(c r > o) — r(c r < o) 
aj = r(c r >o)+r(c r <o)’ 


(3.13) 


is constructed from the partial rates of D decays with Cj greater than or less than zero [155]. 
Correspondingly, the asymmetry for the CP conjugate decay is constructed as 


r (—Cj > o) — r (—Cj < o) 

F(-C r >0)+f(-C r <0)’ 


(3.14) 


where Cj is the triple product for the corresponding CP conjugate particles. 

A non-zero value of the triple product asymmetries aj and Tij themselves is expected due to 
final state interactions (FSI) [155,156]. However, FSI will lead to the same effect for aj and Zij and 
hence their difference is CP asymmetric. This leads to the definition of the CP violating observable 

a ^. p odd = ~( a T -a T ), (3.15) 

where the factor 1/2 simply maintains the definition of an asymmetry. 

Several experimental measurements of triple product asymmetries exist, notably including 
D° —> K K n Ji' and D/. ( —> K®K 1 n n' decays. The most precise results exist for the former 
decay mode, for which the BaBar collaboration measured a T C p° dd = (1.04=5.1 ±4.4) x 10~ 3 [157] 
and the LHCb collaboration measured a T C p° dd = (1.8 ± 2.9 ± 0.4) x 10~ 3 [158]. In particular the 
latest LHCb result shows the impressive robustness of triple product asymmetries against system¬ 
atic uncertainties. The current systematic uncertainty is dominated by the statistical uncertainty of 
a control sample with all other uncertainties being at the level of or less than 10 4 . This robust¬ 
ness is explained in the fact that flavour nrisidentification cannot generate a non-zero asymmetry 
by itself. Furthermore, the resolution on Cj only has to avoid a dilution around Cj ~ 0 and particle 
identification has to avoid confusing the final state products. Since these are not the most challeng¬ 
ing requirements compared to other measurements, triple product asymmetries can be expected to 
remain limited by statistical precision in the foreseeable future and therefore an interesting way of 
searching for CP violation. 


4. Rare charm decays 

Rare decays provide a wide range of interesting measurements. The list of decay modes includes 
flavour-changing neutral currents, radiative, lepton-flavour violation, lepton-number violation, as 
well as baryon-number violation. While a full discussion of rare charm decays would be beyond 
the scope of this paper a few remarks shall be made here. 

4.1 Rare decays and mixing 

There is a direct link between mixing and flavour changing neutral current decays in several exten¬ 
sions of the standard model [76, 159]. These relate AC = 1 annihilation amplitudes, which mediate 
D°—>• / / decays, to AC = 2 mixing amplitudes where the annihilation product, i.e. l~l + , creates 
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c p u 



Figure 9: Re-scattering diagram relating D°—> l l + decays to D°-D° mixing. 


a cu pair from an initial D° meson (see Fig. 9). At tree level, one example is a heavy Z-like boson 
with non-zero flavour-changing couplings. 

The current central values of the D°-mixing parameters, notably x, translate into model- 
dependent limits for rare decays based on common amplitudes. These rare decay limits lie sig¬ 
nificantly below the current experimental limits. D° mixing is well established any upper limit 
from mixing will not change significantly in the future; however, the relevant parameter, x, remains 
to be measured with high precision. Nevertheless, due to the direct correlation of mixing and rare 
decays, any observation above the model-dependent rare decay limits will rule out the correspond¬ 
ing model. The best limit on flavour-changing neutral current decays is the recent LHCb limit on 
the decay D° —y p~ p + of 6.2 x 10~ 9 at 90% confidence level 5 [160]. 

4.2 Rare semi-leptonic decays 

Semi-leptonic decays involving exclusively charged or exclusively neutral leptons proceed via 
flavour-changing neutral current interactions, which are heavily suppressed in the standard model. 
While the decays involving a neutrino anti-neutrino pair are beyond current experimental reach, 
those with a pair of charged leptons can be tested to very high precision. The highest experimental 
precision is achieved in decays involving a p~p + pair by the LHCb collaboration, namely limits 
of 5.5 x 10- 7 for D° —* K^K + p^p + decays [161] and of 7.3 x 10~ 8 (4.1 x 10 7 ) for D + (Df ) —>• 
K + p~p + decays [162]. 

At this level of precision it is necessary to exclude resonant regions that can contribute via al¬ 
lowed decays. Figure 10 shows this at the example of the m(p + p~) spectrum of — y n + p~p + 

candidates. This distribution contains a contribution from —> Xn + decays, withX = (Tj,p, fi), 0), 
followed by an allowed X —> p ji decay. The remaining non-resonant distribution in Fig. 10 is 
entirely due to background processes, but it is in these low -m{p + p~) and high -m(p + p~) regions 
that the search for non-resonant Df-, —> n : p p l has to be carried out. The reported limit on the 
branching fraction is obtained from a limit in these two regions, which has been translated assuming 
a simple phase-space distribution for the non-resonant decay. 

4.3 Lepton-flavour, and lepton and baryon number violating decays 

Among lepton-flavour violating decays the most stringent constraint is a Belle search for D° —>■ 
p e + achieving a limit of 2.6 x 10 7 [163]. Searches for lepton-flavour violating muon or kaon 

5 All further limits are given at 90% confidence level as well. 
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Figure 10: »;(/r + /r ) spectrum of —r n 1 jl p + candidates that pass the final selection and are within 
3ct of either of the mean values of the Gaussian-like peaks describing the D + and l) s signals. Reproduced 
from Ref. [162], 



Figure 11: Feynman diagram of a D A —> n j,L f,i + decay involving a Majorana neutrino. 

decays already provide more constraining limits; however, in scenarious of non-universal couplings 
charm decays, giving access to the up-quark sector, are of great interest. 

The best limit on lepton-number violating charm decays has been placed by the LHCb collab¬ 
oration on the decay D + —> n p p with a limit of 2.2 x 10 x [162]. A particular interest in this 
type of decay lies in its implications on Majorana neutrinos. If the neutrino is its own anti-particle, 
the leptonic D + — y /i 1 V decay can be combined with a virtual anti-neutrino producing a positive 
muon and a negative pion via a boson (see Fig. 11). 

Only the CLEO collaboration has carried out searches for baryon-number violating charm 
decays. Their best limit on the decay D°—> pe~ is 10~ 5 at 90% confidence level [164]. For a more 
complete overview of rare charm decays please refer to Ref. [77]. 

5. Conclusion 

Charm physics provides a vast range of insight into fundamental physics. Clear signals for four- 
quark states including charm quarks have emerged. The width difference of neutral charm mesons 
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is firmly established. 

Nevertheless, a number of open questions remain. It has still not been firmly established 
whether charm mesons can oscillate into their anti-particles or not. Flavour-changing neutral cur¬ 
rent decays remain elusive. And of course the wholy grail of CP violation, well known in the kaon 
and beauty sectors, remains to be found in charm particles despite experimental precision having 
reached levels better than 10 3 . 

With the B factories, CLEO-c and CDF nearing completion of their data analysis, most new 
results will come from LHCb and BESIII. These arc expected to explore very interesting territory 
for charm CP violation. The longer term future will be shaped by the LFICb upgrade as well as 
future e 1 e collider experiments, notably Belle II running at the beauty threshold, and possible 
future charm threshold machines. Charm’s third time has begun to produce its fruits, which may 
well develop into a real charm. 
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